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ABSTRACT 


Convergence zone (CZ) prediction models are developed 
based on acoustic ray tracing theory as applied to linearly 
segmented sound velocity profiles (SVP). The models were 
developed into three calculator programs, two for CZ range 
predictions under different source and receiver depth con- 
ditions and one for CZ gain and transmission loss (TL) pre- 
dictions. The performance of the models as programmed on 
Hewlitt-Packard HP-67 or HP-97 programmable calculators was 
compared to the Fast Asymptotic Coherent Transmission (FACT) 
model which is based on similar but more elaborate theory 
and which is designed for use on large digital computers. 
Agreement of the calculator programs with the FACT model is 
Fairly good when conditions are within the design limitations 


of the programs and environmental conditions are not unusual. 
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፤ 0 ሌሴ ሴ ሌእ ያ እያ ዳጅ ከ ዉስ ው ብክን CZ PREDICTIONS 


፣ ከመው ፡. PROBLEM OF THE AIRBORNE ASW UNIT 

most O the Pacific Ocean and much of the Atlantic 
Ocean and Mediterranean Sea, convergence zone (CZ) conditions 
exist a majority of the time, and they provide passive acous- 
tic sensors with an important means of detecting sounds emit- 
ted from submarine targets. In some areas the CZ regions 
are the most important contact regions considered in planning 
acoustic searches. Obtaining accurate predictions of CZ 
sound propagation is therefore vital to the success of acous- 
tic sensor tactical planners. 

Surrently, the primary source of acoustic predictions for 
U. S. Navy units is the Fleet Numerical Weather Center, Mon- 
terey, California. Propagation loss profiles for four stan- 
dard frequencies and three source and receiver depth 
combinations are normally provided in the ASW Range Predic- 
tion System (ASRAPS) to air ASW units when requested. The 
profiles, showing transmission loss (TL) versus distance, 
are generated on a large digital computer which uses the 
Fast Asymptotic Coherent Transmission (FACT) model. This 
model uses as inputs a linearly segmented sound velocity 
profile (SVP), source and receiver depths, and frequencies 
of interest. The SVP used can be specified by the user or 
can come from information stored at FNWC in the form of 
historical data. This stored data is updated by bathythermo- 


graph (BT) reports through a complex weighting scheme as the 


101 





reports are received. Without going into further detail, it 
can be stated that predictions produced are only as good as 
the data and the computer model used, and only as timely as 
communications allow. 

When entering a search area, a problem often arises con- 
cerning the TL profiles obtained from FNWC. Upon taking a 
BT measurement, the unit often finds the BT profile used to 
generate the acoustic predictions does not agree with the 
Beat BT conditions in the area. If this situation occurs, 
the unit tends to lose confidence in the accuracy of the 
predictions and tactical effectiveness is felt to be dimin- 
ished by lack of good information. The objective of this 
study was, therefore, to investigate what could be done with 
SE or the-art programmable calculators to improve on the 
in situ convergence zone predictions available to air ASW 
Units. 

The reader only interested in the calculator programs 


deyeloped may skip immediately to the appendix. 


RARE LARGE COMPUTERS NECESSARY? 

section 5.6 of Ref. 1, Principles of Underwater Sound, 
by R.J. Urick, the author discusses the relative merits of 
two theoretical approaches to obtaining wave equation solu- 
tions in order to describe the distribution of sound energy 
in space and time. Several references are made to the need 
Memecdigital computers to produce sound propagation descrip- 


tions with either theory. Since those comments were made, 
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however, there has been a revolution in the capabilities of 
meer programmable calculators. Although it is probably true 
that computers are required to produce a complete description 
of sound propagation in the ocean with one program, calcula- 
tors are capable of solving the different modes of propaga- 
tion one at a time with separate programs to obtain a composite 
description. Examples of simple but fairly adequate calculator 
programs for surface duct, bottom bounce, reliable acoustic 
path, and deep sound channel propagation modes are contained 
in Refs. 2 and 3. These references also contain simple models 
memes propagation, but they are based on a mix of ray theory, 
rule-of-thumb, and empirical data. It was felt that a better 
CZ model needed to be developed. 
EENMUDESIRABLE CHARACTERISITCS OF A CALCULATOR PROGRAM FOR 
PREDICTIONS 
i. The program should require a minimum of easily avail- 
able input data. The only information not currently available 
but which would be needed by an airborne ASW unit is an SVP 
from the permanent thermocline to the ocean bottom in the 
ENSE S search area. A chart of this data could easily be in- 
cluded in the environmental package carried aboard the aircraft. 
2. The program should be easy to operate and not require 
the EE uoc have a great deal of insight into the mathe- 
matical model or the internal operations of the calculator. 
3. The output should provide ranges to the inner and 
After edges of all CZ annuli of interest. It should also pre- 
sent the expected TL for all frequencies of interest in each 


annulus. 
T3 





4. The run time of the program should be relatively 
short. This characteristic recognizes that time is impor- 
Bone to an on-station unit. 

5. The program should be based on generally accepted 
acoustic theory. This characteristic is desirable because a 
user would probably have more confidence in such a model 
than one based on empirical data and thus applicable only to 
a specific ocean basin. With empirical models, the user often 
wonders if the area he intends to search corresponds to the 
mean set of conditions used to generate the model or is some- 
how different. 

6. Ideally, the program's performance should agree 
closely with the generally accepted large computer models 


currently in use. 
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II. USING ICAPS AS THE STANDARD FOR COMPARISON 


A. RATIONALE FOR ONLY ONE LOCATION PER OCEAN BASIN 

AS will be demonstrated, the CZ characteristics of the 
three locations studied vary considerably. The deep sound 
channels which produce CZ phenomena are quite different as 
E he ranges from source to CZ annuli. The objective of 
this study was to produce a mathematical model for CZ pre- 
Gictions for use on small programmable calculators. It was 
reasoned that if the model would work for the different con- 
ditions of the three locations studied, it would work for all 
of the variations to be expected within any one of the ocean 
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፣” ድዱ ۴٢٤1 ٨٣ 11+01۱ ርጀጅ ICAPS 

The Integrated Carrier ASW Prediction System (ICAPS) is 
a passive and active acoustic prediction system developed 
WS installation aboard aircraft carrilers and other large 
naval vessels which have digital computers. It contains 
four sets of historical environmental data, one each for the 
North Pacific, North Atlantic, and Indian Oceans, and one 
for the Mediterranean and Black Seas. It also contains sev- 
Pee production programs for predicting naval sonar system 
performance. The FACT model is used in the passive sensor 
predictions. This is the same model used at FNWC for ASRAPS. 
Reference 4 contains a description of the installation and 


Operation of ICAPS in the IBM 360 Computer Center at the 
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Naval Postgraduate School. Reference 5 contains a descrip- 


tron of the mathematics used in the FACT model. 


DESCRIPTION OF HISTORICAL ENVIRONMENTAL DATA FILES 

cures L(a) through 1(1) depict twelve sound velocity 
ENS ues produced by ICAPS from its historical environmental 
data files. Figures l(a) through l(d) show SVP information 
for the months of February, May, August, and November for 
40N 140W in the Pacific Ocean. Figures i(e) through 1(h) 
are the same information for 31N 69W in the Atlantic Ocean, 
and likewise, Figures 1(1) through 1(1) are for 36N 18E in 
the Mediterranean Sea. 

The historical data files used to produce these profiles 
consist of temperature and salinity values for over thirty 
depths, for four seasons of the year, and for many locations 
spaced at one to five degree latitude and longitude inter- 
vals in each ocean area covered. When specific latitude, 
longitude, and date are specified, interpolations are per- 
formed to produce the approximate temperature and salinity 
profiles to be expected at that location and date. This in- 
formation is then converted to an SVP. The output from this 
portion of the system consists of seven columns of values, 
one each for depth in meters and feet, temperature in Cel- 
sius and Fahrenheit, salinity, and sound velocity in meters 
per second and feet per second. The depths associated with 
these quantities begin at ten meter intervals near the ocean 
Surface and gradually increase through 25, 50, 100, 250, 500, 


and 1,000 meter intervals as depth increases. The last line 
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of values is for the ocean bottom depth which was part of 


eae input data. 


E COMPARISON OF DEEP SOUND CHANNEL CHARACTERISTICS 

Figure 2 depicts the deep sound channel (DSC) portion 
of the May SVP for the Pacific, Atlantic and Mediterranean 
coordinates mentioned earlier in a composite graph drawn to 
scale. As can be seen in that figure, DSC characteristics 
of the three areas differ considerably. The vertical ex- 
pent ot the channels varies from 1100 meters in the Mediter- 
ranean to 4200 meters in the Atlantic. The change in sound 
velocity between sonic layer depth (SLD) and DSC axis (point 
of minimum velocity) varies from 15 m/sec in the Mediterra- 
nean to 38 m/sec in the Atlantic. The depth of the DSC axis 
varies from 100 meters in the Mediterranean to 1300 meters 
in the Atlantic. Pacific Ocean values are between the others 
Ball of those characteristics. Sound velocity near the 
surface is much greater in the Atlantic and Mediterranean 
Eu un the Pacific, and sound velocity near the bottom of 
the three basins (not shown in the figure) is about 6 m/sec 
eater in the Atlantic than in the Pacific and about 52 
m/sec greater in the Mediterranean than in the Pacific at 
equal depths. Also note the subsurface sound channel lo- 
cated about 100 to 500 meters below the surface in the At- 


8 ከር profile. 
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NI COMPARISON OF CZ CHARACTERISTICS IN THREE OCEANS 
mm Method Used in Obtaining Data for Comparison 

There were two primary objectives in gathering 
twelve ICAPS runs from each of the three ocean areas. First, 
it was desired to obtain sufficient data to determine which 
CZ characteristics are common to all areas and which charac- 
teristics are peculiar to specific basins. Secondly, it was 
desired to obtain a standard of comparison for any calculator 
program which might be developed. To fulfill the first ob- 
jective, it was decided to keep the input variables the same 
in all areas, varying them one at a time, in order to better 
compare the differences observed in the various runs. For 
each of the three locations, the inputs varied were season 
of the year and source and receiver depth combination.  Re- 
ceiver depths of 60 and 300 feet and source depths of 60 and 
400 feet were used. Each of the ICAPS outputs consisted of 
profiles for four frequencies out to a range of 250 kyds. 

ezıaınally, IE was intended to collect twelve data 
from each profile. These data were to be the range, width, 
McCain, and Transmission Loss for each of the first three 
convergence zone annuli. As it turned out, somewhat less 
data was collected and tabulated. There were several rea- 
Sons for this. First, the February SVP in the Mediterranean 
contained no sound channel and therefore no convergence zones 
existed. Secondly, all of the third CZ data for the Atlan- 
tic was thrown out on the grounds that it was almost always 


the same and that it was inconsistent with information from 


b 





EN ስጎ” ቢሲ two CZ annuli in any particular profile. The 
reason for this occurrence is not know. Finally, it was 
impossible to obtain some of the desired data because of 
the smooth way in which the CZ path blended with other com- 
petitive propagation modes. One could not tell what was 
CZ and what was not in those cases. 
2. CZ Range and Width Analysis 

The transmission loss profiles produced by ICAPS 
are presented in two formats, a table of TL vaiues for each 
kiloyard of range from the source and a graph of the same 
information. Because the TL values are tabulated at kilo- 
yard intervals, it is impossible to be more accurate than 
that interval in determining where a CZ begins and ends. 
Peo, 1t was difficult to be consistent in picking the 
points representing the edges of CZ annuli because of the 
variety of graph shapes, TL levels, and other propagation 
mode interferences. In any event, an attempt was made to 
satisfy one basic criterion in choosing leading and trailing 
edges of the annuli: Do these ranges best represent the 
apparent location of the annulus regardless of the TL levels 
involved?  Admittedly, the ranges picked were often based 
on subjective judgement, and it cannot be stated with com- 
Mete Certainty that only CZ mode propagation contributed 
to the TL peaks judged to be the CZ annuli. 

Table I contains the CZ range and width data that 
Sold be gleaned from the ICAPS profiles. In the table, 


RCZi is the range to the inner edge of the first, second, or 
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۳۶۰۰۰ CZ annulus to the nearest one half nautical mile. CZW 
is the width of the annulus also the nearest one half nauti- 
cal mile. The third column of numbers is the ratio of CZW 
to the range of the outer edge of the annulus (RCZo), ex- 
pressed as a percentage. 

After carefully studying this data, the following 
conclusions were made concerning CZ propagation: 

(1) The range to the first CZ is approximately 
14 to 18 nm at the Mediterranean location, 23 to 27 nm at 
EN 3crfic location, and 33 to 35 nm at the Atlantic location. 

(2) Range to the CZ decreases and annulus 
width increases as source and receiver get deeper in all 
cases. 

(3) The ranges to the second and third annuli 
are approximately whole number multiples of the ranges to the 
inner and outer edges of the first annulus in all cases. 

(4) The range or width of a CZ annulus does 
not appear to have any significant frequency dependence. 

3. CZ Gain and Transmission Loss Analysis 
Convergence zone gain is defined as the difference 
between the transmission loss expected under conditions of 
Spherical propagation and the actual transmission loss ob- 


served. This definition is expressed in Eq. (l). 


Ce= 20, log(r) + atr) = TL (1) 


In this equation, G is the CZ gain, r is the range to the 
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EP Xunnulus, a is the attenuation coefficient associated with 
the frequency of interest, na TL is the actual transmission 
loss observed in the CZ annulus for that frequency. All 
terms in Eq. (1) are in decibels (dB). 

In actual convergence zones, TL (and therefore gain) 
is by no means a constant value. Contributions of several 
possible propagation paths at any one point and the time 
varying nature of sound paths in the ocean cause coherence 
effects to exist. These effects make TL vary in both space 
and time.  Coherence effects are more pronounced at lower 
frequencies (longer wavelengths) where the time varying ef- 
fects are small compared to spatially distributed effects. 
IAPS, the more predictable coherence conditions are in- 
cluded in the mathematical model. 

Since a single TL value was desired for the envi- 
sioned calculator model, an attempt was made to pick the 
average" TL in the ICAPS CZ annuli. As with the range es- 
timates, this called for subjective judgement. Figure 3 
shows a typical ICAPS CZ presentation which has coherence 
effects in evidence. The figure suggests how an "average" 
TL was chosen as best representing that E TWO levels 
were chosen (labeled high and low in the figure) which bracket 
Piewmajority of the TL points within the annulus. The ap- 
proximate midpoint between those levels was then picked as 
meme TL for that CZ. ` 

AS an extra point of interest, the high and low TI 


levels were studied. It was noted that ICAPS predicts TL 
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variations from about Y 104B around the "average" level. If 
this is truly representative of CZ coherence effects, an ASW 
unit armed only with an estimate of the "average" TL in a 
certain CZ annulus should expect to see variations of about 
mee magnitude around the estimate in hand. 

Using TL levels estimated by the procedure described 
above, and employing Eq. (1), CZ gain values predicted by 
ICAPS were obtained and tabulated. The values produced are 
contained in Tables II and III. Table II shows all of the 
Gee from the Pacific location. Table III contains only 300 
Hz data from the Atlantic and Mediterranean locations. (50, 
850, and 1700 Hz data were omitted from Table III because it 
became obvious during data collection that G is not frequency 
dependent.) 

Again after careful study, the following conclusions 
were drawn concerning CZ gain: 

(1) CZ gain values range from about eight to 
twenty dB in all three areas observed. 

(2) CZ gain is the same value for first, 
second, and third CZ in any given case. 

(3) In general, CZ gain is independent of 
frequency. An exception to this conclusion is that at low 
frequency (below 300Hz), especially when source or receiver 
or both are above the SLD and/or near the surface, there is 
apparently somewhat less gain than evident for higher fre- 
quencies. This difference is probably due to stronger dif- 


fraction of the longer wavelengths. 
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(4) CZ-.gain seems to be highest when source 


anda receiver are at or near the same depth. 
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Month REV eat Freq 02 Gain dB) 


een Tr 15ቲ eg na 02. 3ሯሷ CZ 
50 7 7 6 
300 14 13 12 
00/60 350 14 6 E 
1700 1 E 12 
50 9 9 9 
FEB 300 12 12 iL 
246 60/ 400 850 11 11 12 
1700 11 i 12 
50 12 11 10 
300 12 12 11 
300/400 850 10 il 12 
1700 12 10 ii 
50 16 1 de 
300 ic 17 16 
60/60 350 Te 17 16 
1700 We 17 l7 
50 10 10 12 
MAY 300 10 5 js 
33 60/400 850 el 12 15 
1700 9 il 12 
50 14 14 14 
300 E 13 14 
300/400 850 e 14 iu 
1700 12 EI = 


Maple II. ICAPS CZ Gain Data Observed at the Pacific Ocean 
Location. (Page 1 of 2) 
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Month Rewer Pg t Freq CZ Gain (dB) 





SLD(ft) (ft) (Hz) Ist cz Ind CZ 3rd CZ 
50 11 11 13 
300 [4 13 12 
60/60 850 e 18 20 
1700 16 Ts 18 
50 7 7 9 
AUG 300 9 12 13 
0 60/400 350 12 14 1٠ 
1700 12 13 6 
50 14 14 14 
300/400 300 11 14 12 
850 e 12 E 
1700 12 LE 11 
50 3 7 6 
300 9 12 12 
60/60 350 15 18 P 
1700 14 "s 15 
50 9 9 9 
NOV 300 12 Us i 
98 60/400 850 11 11 1 
1700 11 11 11 
50 11 16 1 
300 12 11 m 
300/400 850 9 3 12 
1700 10 9 7 





Table II. ICAPS CZ Gain Data Observed at the Pacific Ocean 
Location. (Page 2 of 2) 
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Ocean Month ከ IgE O AMO da a 





SLD(ft) (Et) lst CZ ond CZ 3rd CZ 
60/60 17 14 
-= 60/400 13 13 
300/400 (Es 18 
60/60 19 14 
== 60/400 16 12 
5 300/400 13 17 
» 
E 
E ٣ 60/60 15 17 
E = 60/400 10 TE 
300/400 14 6 
60/90 15 11. 
-፦- 60/400 12 10 
300/400 13 12 
60/60 14 Ls d 
MAY 
m 60/400 12 10 10 
300/400 15 18 20 
52 
S 60/60 17 19 1.7 
AUG 
< I 60/400 E 9 10 
z 300/400 
t 
E 
2 60/60 11 e ገር 
NOV 
= UE 60/400 10 m 12 
300/400 i 14 E 


Table III. ICAPS 300Hz CZ Gain Data Observed at Atlantic 
and Mediterranean Locations. 
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III. CZ RAY THEORY ANALYSIS AND MODEL DEVELOPMENT 


A. CZ RANGE AND WIDTH 

It was decided to use ray tracing as the method for de- 
termining CZ range and width because of the simplicity of 
the mathematics involved and because of the intuitive appeal 
Reg rays depicting the propagation of sound. The alter- 
native approach, that of normal mode theory, was rejected on 
the grounds that it would be much more complicated, requiring 
capabilities far beyond those available in the calculators 
at hand. 

Figure 4 shows four sound Hays Or Particular interest in 
CZ propagation. The order of these rays is described for the 
"typical" case in the following discussion. An "atypical" 
case will be mentioned later. 

Ray #1 departs the SLD at zero degree depression angle. 
Įt reaches its greatest depth at the bottom of the DSC and 
returns to the SLD at some particular range and at horizontal 
incidence. The horizontal range from SLD to SLD is termed 
cycle distance. The cycle distance for this ray is desig- 
nated Co: 

[p Zus che next ray of interest found as the departure 
angle from the SLD is increased downward. This ray passes 
down through and below the bottom of the DSC before turning 
back upward. It returns to the SLD at the shortest range 


from the starting point of any ray within the bundle of rays 
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undergoing CZ refraction. This range is designated rnin 


The angle of departure for this ray is designated ee 
Each ray between £l and #2 crosses all of the previous (les- 
ser departure angle) rays on its way up from its lowest 
depth. 

As departure angle from the SLD is further increased, 
the next sound ray of interest, $3, is located. This ray 


has a cycle distance equal to r Its maximum depth is 


o: 
greater than that for ray #2. It departs from and arrives 


back at the SLD at an angle designated 9 Rays between 


rswp' 
#2 and #3 do not cross each other, but they do cross the 
earlier rays on their way back up to the SLD. 

In the CZ annulus, the rays between #1 and #2 sweep in- 
ward toward the source as departure angle increases. After 
ray #2 they sweep out away from the source as angle increases 
further. For this reason, the region formed by rays between 
#1 and #3 is called the reswept zone. 

Finally, as angle of departure from the SLD is increased 
to maximum angle for CZ propagation, we observe ray #4. This 
ray turns upward at a depth equal to the water column depth 
at that location. It returns to the SLD at the greatest dis- 
tance of all CZ refracted rays. Rays departing the SLD at 
angles greater than that for ray #4 would be reflected off 
the bottom and are not of significance for the CZ propagation 
path. 

As mentioned earlier, this progression of rays exists 


in a "typical" CZ situation. If, however, the ocean bottom 
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were more shallow, cycle distance for ray #4 would be re- 
duced. If the bottom were shallow enough, ray #4's cycle 
distance would be less than Ca: In that case, the reswept 


zone would be reduced to the region between rays #2 and #4. 


This situation is called the "atypical" case. 


B. RANGE AND WIDTH MODEL 

In the calculator programs developed, provision is made 
femme ntering and storing a five point sound velocity profile 
which defines the DSC only. The first depth and velocity 
pair entered (0317) ) equate to the appropriate values found 
at the SLD. The fifth depth entered (Dz) is the depth at 
the bottom of the DSC where sound velocity is equal to that 
at the SLD. The other three depth/velocity pairs must be 
picked subjectively from a graph of the SVP of interest. If 
a mixed layer exists, the gradient in that layer is taken to 
be 0.02 sect (a purely pressure induced gradient to two 
Place accuracy). The program calculates the four layer gra- 
memes within the DSC profile entered, and uses the fourth 
(deepest) layer gradient in ray calculations that occur below 
the DSC. It would have been desirable to allow several more 
points in the SVP, but calculator data storage capacity and 
program step limitations preclude more than five depth/veloc- 
NE pairs. 

The overall scheme used to predict CZ range and width is 


tO trace a series of rays starting at the SLD with a zero de- 


pression angle ray. That first ray yields rg which is stored. 
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Then an iterative process is begun in which the angle is in- 

cremented and each succeeding cycle distance determined is 

Reese to the previous one untilr . and 9. .. are found. 
min rmin 


(8 is stored for use in the CZ gain and TL program, to 


rmin 
be discussed later.) Corrections are then made to د1‎ and 
Lo to account for surface duct effects (if any) and source 
and receiver depth separation from the SLD. Range to the 
inner edge of the CZ is Ee plus corrections, and range to 


the outer edge is r, plus corrections. 


0 
ihe first attempt to produce a calculator program, 

the cycle distance of ray #4 (the ray just grazing the bot- 

tom) was compared to ማሽ The greater of the two was picked 

as the basic distance for determining range to the outer 

EuSGOOf a CZ. Later on, this portion of the program had to 

be deleted to save program steps. The final programs de- 

veloped ignore bottom depth and do not include rays outside 

the reswept zone in determining annular width. This 1s prob- 

ably a shortcoming of the programs but the seriousness of 

the errors it causes will not be known without further study. 
A commonly applied rule-of-thumb states there must be 

a minimum 300 fathoms of depth excess (water column below 

the DSC) in order to have "reliable" CZ conditions. It was 

observed that a fully developed reswept zone existed in 

every case in the locations studied, and separate calcula- 

tions showed that somewhat less than 300 fathoms depth ex- 


cess was required to complete the zone. Therefore, a program 


user should consider the 300 fathom rule-of-thumb before 





running the range and width program. With less than 300 
fathoms depth excess, the possibility exists for an "atypical" 
CZ propagation situation where the reswept zone is reduced in 
width due to bottom ray limiting. 

Another program shortcoming involves an assumption that 
both source and receiver would be more shallow than the sec- 


EN DSC SVP point chosen (depth D,). In other words, the 


2 
programs were designed to allow for source/receiver depths 
within the mixed layer or the first isogradient layer below 
the SLD. After the five point SVP is entered and the grad- 
ients computed, source and receiver depths are entered and 
converted to velocities. The programs determine these veloc- 
ities (Ce and Cp? by subtracting an appropriate amount from 
the velocity at the SLD. The amount subtracted is determined 
by depth separation from the SLD and by the gradient in elther 
the ML or the first layer below the SLD. If source or re- 


ceiver depth 1s greater than D sound velocity should be de- 


SE 


termined by correcting C, (the velocity at 22) and by using 


2 
وو‎ (the second layer gradient). Since this is not done, ve- 
locities for source/receiver depths below D^ will belin error 
(usually too low). Source and receiver velocity errors are 
Carried over into Ar. and Ar, range corrections. If the 
velocities are too low, the range corrections will be too 
Marge. This is normally a rather insignificant source of 
total range error, however, since Ar. and Ar, errors will be 


a Small fraction of the magnitude of those terms and because 


the range correction terms are small to begin with. 





The mathematics of ray tracing in isogradient layers is 
quite straightforward. By Snell's Law, 


a z 2 Bee ስሮ ከ 
cos 94 COS SÉ for each ray) 


(2) 


the angle of a ray departing a layer can be determined from 


Enola Of entry into that layer. In Eq. (2), C, is the 


1 


sound velocity where the sound ray enters a layer, 91 is the 


angle of entry, C, is the sound velocity where the ray de- 


2 


parts the layer, and 9, is the angle of departure. 


ves Eravel 11 Circular afcSs within constant gradient 
layers, and the radius of curvature is: 


a 


g, cos 91 


(3) 


L 1[ 


the layer (in this case, layer ۰ 


C, and 9, are as defined above and 21 is the gradient within 


Finally, the horizontal distance traveled bv a ray 

while traversing a layer is: 
Ar = DIG 25 SAA 97) (4) 

Figure 5 demonstrates an example application of Eqs. 
2 through 4. It should be noted that absolute value signs 
are used in Eqs. 3 and 4 because the gradient in Eq. 3 and 
the difference of sines in Eq. 4 may be positive or negative, 
while R and Ar are always positive. 

The programs use these equations to compute the horizon- 
tal range increments each ray accumulates within the four 


layers, doubles each term (to account for the downward and 
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Sound Velocity ሮሙ Range 





Figure 5. Horizontal Distance Traveled Within 
an Isogradient Layer. 
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upward passes through each layer), and then sums the terms 
to obtain cycle distances. The fourth layer requires a 
Slightly different treatment because the rays become horizon- 
tal and then turn back upward within that layer. Essentially 
the same formulas are used, however. The equations are also 
used to compute the range correction terms. 

In considering the various possible ray paths between 
source and receiver, it was decided there were four basic 
ENEDations which could occur: 


1) No mixed layer, both source and receiver 
below the SLD.  (Deep/Deep) 


2) Mixed layer present, both source and receiver 
below the SLD. (Deep/Deep/ML) 


3) Mixed layer present, both source and receiver 
within the layer. (Shal/Shal) 

4) Mixed layer present, source or receiver above 
the SLD, the other below. (Crosslayer) 


The only difference between the first two cases is the mixed 
layer effect in case 2, That effect causes a widening of 
annuli due to spreading of sound rays as they travel up to 
the surface and back down to the SLD within the layer. The 
mixed layer effect is also included in the third and fourth 
cases above. 

It was originally intended to include all four cases in 
one range prediction program. Again due to calculator limi- 
tations, it was necessary to use two programs to cover the 
four possibilities. The first range program (labeled Deep/ 
Deep) is for cases 1) and 2) above when both source and re- 


eeiver are below the SLD whether or not an ML exists. The 


25 








second range program (labeled Shal/Shal or Crosslayer) is 
for use in cases 3) and 4) above when source or receiver or 
both are above the SLD. 

Formulas used to determine range to inner edge of the 
first CZ (RCZi) and range to the outer edge of the first CZ 


REZO) follow: 


RCZi - ۰ e Arg - Ara Deep/Deep 
د‎ + Ar. + Alo Shal/Shal 
ERE ር Ar. 5 Ar. Crosslayer 
RCZo = r, + f, e, uu A گ تت د‎ 
^ rg + 2 Arg - Arg - Ar. Shal/Shal 
= ro t 2 Ary E Ar. T Ar. Crosslayer 
In these equations, Dem and Co have been previously defined, 


Ar, and Ar, are the respective horizontal range corrections 
which account for source and receiver depth separation from 
Ee SLD, and 2 AYg is the correction toy mixed layer effect. 
eres 6(a) through 6(f) (not to scale) depict the RCZi and 
RCZO formulas in graphic form. Ranges to second and subse- 
quent CZ annuli are taken to be integer multiples of the 
ranges produced. 

Two items of interest, both evident in Figs. 6(a) - 6(f), 
More mentioning at this point. First, acoustical reci- 
procity is envoked and the more shallow of source and re- 


ceiver is always treated as "source" of the sound rays within 


pres calculator programs. Secondly, only those sound rays 
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Figure 6(a). RCZ, for Deep/Deep Case. 
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Figure 6(c). RCZ, for Crosslayer Case. 
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which experience no more than one ocean surface reflection 
between source and receiver are considered in this model. 
Both of these conventions are commonly applied to ray trac- 
ing models. Although they theoretically have’ little or no 
effect on model results, they greatly simplify the work of 


programming a ray tracing model. 


۰۰۰.٠7م‎ GAIN AND TRANSMISSION LOSS MODEL 
In general, transmission loss is defined as ten times 
the logarithm of the ratio of sound intensities measured at 


One meter from a source and at range r from that source. 


= 
EE <> IO iLexe; Ee 
ገር 
Intensity has units of power per unit area. The change in 


intensity between one meter and range r is due to geometric 
spreading of the power over a different amount of area and 
due to attenuation of some of the power through absorption, 
scattering, diffusion, etc. 

In ray tracing theory, it is assumed there is no sound 
power transfer across sound rays. Therefore, the power flow- 
ing from a source between two sound rays remains between 
those rays and travels out in a direction parallel to the 
ray paths. Determining the portion of transmission loss due 
to geometric spreading I under this assumption reduces to 
mending ten times the logarithm of the ratio of areas (at 
range r and at one meter) penetrated by the power between 


the two rays perpendicular to the direction of travel. 
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TM NDS Tog —— 
g 72 ጺ) 


A mathematical development of this technique is contained on 
pages 119-121 of Ref. Il. 

Figure 7 shows how this method was adapted for use in 
the CZ Gain and Transmission Loss portion of the calculator 
model developed. The area (Ay) at one meter from the source 
is the product of area height and area circumference. The 
sound rays bounding the area above and below are the minimum 
and maximum departure angle rays which produce the reswept 
zone in the CZ annulus. The angular spread of those rays 
(A48) in radian units times the sphere radius (l meter) is 
the area height. Cosine of the average angle of departure 


of the rays (0,) times the sphere radius times 27 is circum- 


1? 


ference of the area. Therefore: 


E : 2 
Ar = 2m ÁO cos 91 (m) 


Tn the CZ, the area (A,) over which the same power is 
Peeseeibuted 1S also found by a product of area circumference 
and area width. Circumference is 27 times range to the CZ 
(RCZi). Width of the area perpendicular to the sound rays 
mm hie product of CZ annulus width (CZW) and the sine of the 
average angle of arrival of the rays at the receiver depth 


(9 Therefore: 


2) * 


2 


ወ memes CZW Sin 9, onm 


2 


and the geometric TL expression becomes: 
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RCZi CZW sin 9 


2 
A8 cos 9 


ieee 19 
2 1 


Substituting this expression back into Eq. l, which is 
the definition of CZ gain, and reducing to simplest form 
yields the algorithm used to determine G in the calculator 
model: 


REZ A0 COS 91 (5) 


CZW sin 95 


G = 10 log 


To implement this algorithm, the program has only to 
determine the angular terms since RCZi and CZW are available 
from the range program results. After one of the range pro- 
grams has been run, the user loads the G and TL program into 
calculator memory without altering the contents of the data 
storage registers left from the range program. Then the 
iterative ray tracing process begun in the range program is 
continued in the gain program. The angle of departure of 


sound rays from the SLD is incremented beyond 9 (Len 


Ena 


storage from the range program) and cycle distances produced 


are checked for approximate equality with r In this way, 


Q^ 
the ray which completes the reswept zone is found, and its 


angle of departure from the SLD is TS Then Eos and 


zero degrees (the angle for the ray producing cycle distance 


i are converted to angles of departure from the source 


0) 
depth, Dan and 90 respectively, and angles of arrival at 


the receiver depth, Den and 9 50 respectively, using Snell's 
law. The angular terms in Eq. 5 are then computed using 


the following formulas: 
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A0 26 ET 


SR SO Deep "source" 
= 0 چو‎ 0 " " 
SR SO Shal "source 
Ee "SR E "SO Deep "source" 
1. ረ 
E ISR 0 Iso Shal "source" 
2 
T "an" Tas 
2 2 


Recall that "source" in the model refers to the more shallow 
Of source and receiver. Therefore, the deep "source" forms 
of these formulas are used only after using the Deep/Deep 
Busse program. In all other cases, the "source" is consid- 
ered to be shallow. Figure 8 depicts these angular relation- 
Pees or the various depth conditions. 

It should be pointed out there are two inherent errors 
in the angular quantities determined. First, the possible 
source and receiver sound velocity errors mentioned earlier 
could cause the gain algorithm angles to be slightly off. 


This would only occur if depths greater than D, were entered 


2 
for source or receiver or both. Secondly, s is found 
for the ray which has cycle distance equal to Co at the SLD. 


Since the actual CZ ray bundle departs from the source depth 
(vice SLD) and arrives at the receiver depth (vice SLD), the 
ray which completes the reswept zone will probably be differ- 
ent than the ray used and it will have a slightly different 


angle crossing the SLD. These angular errors will cause the 


6.3 
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Figure 8. Determining angular terms for 


Ca ain algorithm. 
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greatest CZ gain error in the sin SÉ term of the algorithm. 
Since sine is directly proportional to angle at small angles, 


ከ የን" ርቶ ሟ factor of two in 9 (a quite possible event) 


2 
could cause a gain error of approximately 3 dB. 
Once CZ gain is computed and stored, the sound frequency 


Of interest is entered, and the attenuation coefficient is 


ሬክ  ቢር0ዐ| using Thorpe's equation (p. 102, Ref. 1): 


2 2 
ው (0.001094) + —À (dB/m) (6) 
l + £ 4100 + £ 


EM 6, f Is in kHz, and the constant in front of the ex- 
pression converts attenuation coefficient from dB/kyd to 
dB/m. The program user enters frequency in Hz, and the pro- 
gram performs the conversion to kHz. 

th 


Finally, the transmission loss in the n CZ annulus 


for the frequency of interest (TL, ) is determined by: 


)7( ته خو( دت د TL, SS A‏ 


E 6۰-27: 101:07 


In this equation, the subscript n denotes the n 
the range to which is n times RCZi. 

After a range program is run, and after the gain portion 
of the G and TL, program has been completed, TL. values for 
a variety of frequencies and CZ annuli may be rapidly ob- 
tained for the SVP, source depth, and receiver depth condi- 
rons entered. If, however, a different set of source/receiver 
depth conditions are also of interest, the entire procedure 


Beginning with the appropriate range program must be per- 


formed again. 


65 





ዘ መይ CULATOR PREDICTIONS COMPARED TO ICAPS 
۳۰۰٠٠۰7٦۰ tne SVP Points for the Program 

The five point SVP limitation of the ray tracing 
procedure is a rather serious handicap in many situations. 
Actual sound velocity profiles not only are curvilinear in 
overall shape but also have many small scale features and 
they are time varying functions. Approximating these curves 
with only four straight line segments presents a difficult 
challenge. 

In general, matching the gradients, sound velocities, 
and associated depths are all important in choosing SVP 
points. The greatest potential for causing large range pre- 
diction errors occurs when the SVP contains an extensive 
near surface layer with very slight velocity gradient. The 
horizontal distance traveled by a shallow depression angle 
ray within such a layer varies considerably with small 
Changes in the gradient or layer thickness. Under such con- 
ditions, then, it is extremely important to match those 
characteristics as closely as possible. 

Another important item to carefully match is the 
sound velocity at the DSC axis. This velocity determines 
the maximum angle of depression for each ray prior to com- 
mencing upward refraction. The horizontal distance traveled 
by a ray below the axis is highly dependent on that angle. 

Table IV contains the five point sound velocity 
profiles picked by the author for use in comparing the pro- 


gram performance to ICAPS predictions.  Depths in the table 
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are in meters, and sound velocities are in meters per second. 
The reader may want to plot these points on the graphs of 
Figs. l(a) through 1(1) so he may see how the four isogradi- 
ent layers picked match the ICAPS profiles. It should be 
pointed out that only the initial selection of SVP points 
was used in the subsequent comparisons of calculator model 
results to ICAPS predictions. Since an ASW aircrewman using 
the programs in attempting an in situ prediction of acoustic 
conditions would not be able to judge whether SVP point ad- 
justments would improve or degrade prediction accuracy, it 
was felt that comparing resuts of the initial SVP point se- 
lection with ICAPS would be more meaningful to the objective 
of developing the calculator programs. 

Coneasing calculator model predictions to ICAPS pre- 
dictions in a definitive statistical manner was not done. 
The main reason for this was alluded to in the preceding para- 
٣۰٠۰٠۰٦٠۰ Since the SVP points entered in the calculator pro- 
gram must be picked subjectively by the person using the 
program and since it is unlikely different people would pick 
the exact same points off any given SVP, it is clear that 
Calculator results can be expected to vary from operator to 
operator. 

2. CZ Range and Annulus Width Comparisons 

The calculator range programs produce one value each 
Bor RCZ] and RCZo for any given SVP, source depth, and re- 
ceiver depth situation. Under the same set of conditions, 


ICAPS yields four sets of RCZi and CZW values, one set for 
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E. of the four frequencies entered. In order to compare 
the calculator performance to ICAPS it was first necessary 
to reduce the ICAPS predictions to one value each for ጀርመ 
and CZW for each SVP/source/receiver condition. This was 
done by simple averaging to eliminate the frequency variable 
from the ICAPS range and width predictions. 

Figures 9(a) - 9(c) display range and width compari- 
eons in graphical form for the Pacific, Atlantic, and Medi- 
terranean locations respectively. In each figure the double 
barred lines represent the ICAPS first CZ annuli predictions 
(averaged over frequency), and the single barred lines repre- 
sent the calculator predictions. Numerical values for inner 
ከደ ከ ር: first CZ ranges may be obtained from the scales at 
the tops of the figures. 

In all, there were 32 cases where these graphical 
comparisons could be made. The following comments pertain 
to tnose comparisons: 

a) In 30 of the 32 cases the calculator annuli 
overlap at least a portion of the ICAPS annuli. 

b) In 14 of the 32 cases the calculator annuli 
are completely contained within the limits of the ICAPS 
annuli. 

C) In all 32 cases RCZi ranges predicted by 
the calculator were greater than those predicted by ICAPS. 
Ehe 12 Pacific cases, the calculator RCZi values were ap- 
proximately 1.7 nm greater than ICAPS on the average. In 


the 12 Atlantic cases, the average difference was approximately 
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Figure 9(a). Comparisons of ICAPS and Calculator CZ 
Annuli Predictions. Pacific Ocean. 
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Atlantic Ocean. 
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Figure 9(c). Comparisons of ICAPS and Calculator CZ 
Annuli Predictions. Mediterranean Sea. 
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2.8 nm. And in the eight Mediterranean cases, 1.5 nm was 
the mean difference. 

d) In 27 of the 32 cases the CZ width predic- 
tions from the calculator were more narrow than the ICAPS 
predicted widths. Three of the five cases where calculator 
CZW exceeded ICAPS CZW were from the February SVP in the At- 
lantic location. That SVP contained a very deep (500 meter), 
nearly isovelocity layer near the surface. In such a pro- 
file, CZ refraction produces ray paths that are spread over 
a very wide (in this case 16-20 nm) annulus. Only the rays 
which return to the SLD within the first few nm at the inner 
edge of that annulus experience sufficient convergence to 
produce detectable CZ gain, however. Going from inner to 
outer edge of such an annulus the CZ refracted rays rapidly 
fan out experiencing progressively less convergence and pro- 
ducing progressively less CZ gain. Additionally, if the 
bottom grazing ray were considered, it would be seen to limit 
the reswept region of this type annulus to something far less 
Em Ehat indicated. Since the calculator model fails to ac- 
count for either of these factors, it fails rather dramat- 
MERA to produce a "practical" CZ annular width from this 
EN DES, 

In summary, the calculator model produces CZ annuli 
that roughly agree with those produced by ICAPS in all three 
Ocean basins considered. Calculator RCZi ranges are 5-10$ 
greater on the average than the ICAPS values. Calculator 


CZW predictions (excluding the Atlantic February SVP) are 
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40-50% narrower than ICAPS widths on the average. And, the 
Atlantic February case indicates there is at least one SVP 
type in which the calculator model fails to produce even 
marginally acceptable results for CZW. 

3. CZ Gain Comparisons 

As with CZ range and width comparisons, it was nec- 
essary to average the ICAPS gain data with respect to fre- 
quency before calculator gain predictions could be compared. 
The estimated ICAPS gain values in Tables II and III were 
Buusereduced to one number for each SVP, source, and receiver 
condition. Table V, CZ Gain Prediction Comparisons, contains 
numbers that represent the difference between calculator gain 
predictions and the averaged ICAPS values. Minus signs in 
the table indicate those cases where calculator gain was less 
tnan the ICAPS value. 

As with the range and width comparisons, the worst 
agreement occurred in the Atlantic winter SVP case. Since 
CZW is a term in the gain algorithm, the extremely wide an- 
nuli predicted by the calculator caused gain values to be 
far too low for the three source/receiver conditions assoc- 
lated with that SVP. 

Excluding the Atlantic winter SVP case, the follow- 
ing comments can be made concerning the other 29 CZ gain 
comparisons: 

a) Calculator gain values ranged from 7.3 dB 


lower to 5 dB higher than the averaged ICAPS values. 
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ESE v. CZ Gain Prediction Comparisons. 
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E ST Ele 29 comparisons, calculator 
values were within 3 dB of ICAPS. 

c) In nine of the 29 comparisons calculator 
values were within one dB of ICAPS. 

d) On the average, calculator gain values were 
approximately one dB less than ICAPS. This result is incon- 
Sistent with calculator CZW results in light of the gain mod- 
el used. Since calculator CZW values averaged only slightly 
more than half the ICAPS widths, 1t would have been more con- 
Sistent 1f calculator gain values turned out two to three dB 
higher than ICAPS (acoustic power being spread over less area 
in the CZ annuli, other things being equal). Perhaps an ex- 
planation for this apparent discrepancy is that the calcula- 
tor model does not consider the contribution of surface 
reflected energy adding to the energy from upward traveling 
sound rays at the receiver depth. In an actual CZ annulus 
the downward traveling, surface reflected energy adds approx- 
imately three dB to the CZ gain over much of the annulus 
width. Apparently, the FACT model in the ICAPS system in- 
Eves this consideration. It is also apparent that neglect- 
ing surface reflected energy in the calculator gain model has 
the effect of canceling errors that should result from CZW 
values being too narrow. 

In summary of the gain results, it can be said that 
the ray tracing technique used in the calculator model worked 


reasonably well. Since three fourths of the comparison cases 
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resulted in gain values within three dB of the estimated 
ICAPS figures, TL values from the calculator displayed the 


same close agreement. 
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۶9 CONCLUSIONS 


IMITATIONS OF THE MODEL DEVELOPED 

The HP-67/97 calculators used in programming the CZ pre- 
diction model were stretched to their limits in both data 
storage and program step capacity. Although not known for 
certain, the author feels significantly more accurate results 
would be possible from a calculator with only moderately 
larger storage capacity. 

The data storage limitation which allowed only five SVP 
points to be entered is quite restrictive and no doubt plays 
a large role in the CZ range and width inaccuracies obtained. 

Program step capacity forced several short cuts to be 
taken which again would not have been necessary with a mod- 
erately larger program memory. Two separate range and width 
programs were required due to insufficient program space to 
incorporate tests for different source and receiver depth 
cases. Also, source and receiver depths are strictly allowed 
only within the upper two SVP isogradient layers because pro- 
gram space was not available to check for the correct layer 
if all depths were allowed. Additionally, and perhaps the 
greatest source of CZW errors observed, program step limita- 
tion prevented incorporating a method of considering the 
bottom limited CZ sound ray in determining the range to the 
outer edge of a CZ annulus. The program developed ignores 


the bottom entirely and considers only the reswept zone in 


78 





predicting annular width. Since calculator CZW results were 
Pemsiderably shorter than those indicated by ICAPS, it is 
assumed the discrepancy is due to not considering CZ rays 
beyond the reswept zone. The first priority in making im- 
provements to the calculator model, should a larger capacity 
machine be implemented, would be incorporating a better meth- 
Bomtor Selecting the ray which defines the outer limit of 


mer EZ annulus. 


IIS EFULNESS OF THE MODEL DEVELOPED 

The degree of success in producing a useable CZ prediction 
model for handheld calculators must be determined by consid- 
Meme une objectives set forth in the first section of this 
Peay, The Central idea was to ascertain if a calculator 
Enc would improve on ASRAPS CZ prediction accuracy in the 
case where BT conditions determined in situ differed from 
those used to generate the ASRAPS TL profiles.  Inherent in 
this objective is the assumption that ASRAPS TL profiles 
generated primarily from climatological data would be in 
SS CU CO lack of input data accuracy. Also inherently 
med 15 that given identical input data the calculator 
model would produce less accurate results than the digital 
computer model(FACT) used in ASRAPS(and ICAPS) due to obvious 
differences in data and program capacities. The real question 
then is a trade-off comparison: Will the basically less accu- 
rate calculator model produce better CZ predictions with actual 
environmental data than the more sophisticated digital com- 


puter model which had only climatological input data? 
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Before addressing the answer to this question, character- 
istics of the calculator model developed will be compared to 
the list of six desirable characteristics described in sec- 
morn T. 

e easily available input data. 

The data required are an SVP, assumed source 
depth, hydrophone depth, and frequency of interest. The 
only portion of this information not presently available to 
ASW aircrews is that part of the SVP below the 1,000 ft 
depth limit of the AN/SSQ-36 bathythermograph buoy. SVP data 
from the surface to 1,000 ft (the area where seasonal and 
diurnal variations predominantly occur) is easily obtained 
Beem che BT buoy information. 

2. Ease of program operation. 

an OR. ۹3211113 11 HP=6//97 calculator use 

could operate this program without additional training. 
ይ ۲006 038 ۔‎ 

The program provides CZ annulus range and width 
as well as TL values for all frequencies of interest in all 
Mull of interest. 

T ES Hope ۲ 1۷۹۷717167 

TO run a complete program requires approximately 
10 minutes once SVP data is obtained. Deploying a BT buoy 
and converting the temperature trace to an SVP would take an 


additional 10-15 minutes. 
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ሰጠን ታያ entirely on acoustic theory. 

The program uses only ray tracing techniques in 
producing its output terms. 

6. Agreement with large computer models. 

Calculator CZ ranges obtained averaged 5-10% 
greater than ranges obtained from ICAPS. CZ width results 
averaged only 40-50% of those obtained from ICAPS. And, 
there was one SVP case studied (winter, Atlantic) in which 
the calculator CZW results were very different from ICAPS. 
That SVP case was considered a failure of the calculator 
model, and it must be conceded the model does not work for 
AAN Situations. Excluding the obvious CZW failure SVP 
case, TL values from the calculator averaged about one dB 
lower than ICAPS with extreme deviations observed ranging 
from -7.3dB to +5dB around the ICAPS values. Additionally, 
Calculator results can be expected to vary from operator to 
operator since SVP points must be picked subjectively from 
mM SVP graph. 

Returning to the main objective of the study, the author 
feels that only half of the trade-off question has been 
answered. An easily operated, purely theoretical model was 
developed which works for most (but not all) CZ producing 
SVP conditions. And a measure of its accuracy compared to 
the sophisticated FACT computer model was obtained. Yet to 
be answered is how inaccurate ASRAPS CZ predictions are when 
boserved BT conditions differ from climatological conditions. 


This portion of the question is very difficult to answer and 


5.1. 








in fact would be a very large study in itself. Generally, 
inaccuracies must range over a scale from insignificant to 
considerable as environmental deviations range from slight 
to great. The most likely variables affecting degree of in- 
accuracy are surface water temperature, mixed layer depth, 
and in layer andbelow layer gradients. The effects of vary- 
ing these or other possible factors one at a time or in 
various combinations on CZ range, width, and gain must be 
known before the entire question can be answered. Further, 
definite magnitudes of environmental factor deviation must 
be determined so that a person can judge when ASRAPS inaccu- 
racies are likely to be greater than the calculator model 
inaccuracies. Until these points are answered it would be 
inappropriate to recommend use of the calculator model as a 


routine method of updating ASRAPS CZ predictions in situ. 
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APPENDIX 


HP-67/97 Calculator Programs for Convergence 


zone Range, Width and Transmission Loss Predictions 


Steps required to use the programs: 


l. Deploy a bathythermography buoy in the operating area 
of interest. 


COn Vert the BT buoy information to a sound velocity 
profile of the upper 1,000 ft of the ocean area. 


3. Combine the upper SVP data with a graph of climato- 
logical SVP data which depicts sound velocity conditions 
below the 1,000 ft level. 


4. Pick five points from the combined SVP graph which 
best represent the deep sound channel portion of the SVP. 
ከመ ist point should be at the sonic layer depth, the fifth 
point at the bottom of the DSC where sound velocity equals 
mat at the SLD, and the other three points at points on the 
graph such that when straight lines are drawn to connect the 
five points they create a linearly segmented SVP which matches 
the actual SVP as closely as possible. 


5. Pick the appropriate CZ Range and Width program to be 
used as follows: If both source and receiver are below the 
SLD use the Deep/Deep program. If source or receiver or both 
are above the SLD use the Shal/Shal or Crosslayer program. 


6. Load and run the appropriate Range and Width program 
Peeenaing to the accompanying instructions. 


Ss ina the calculator power on and data storage reg- 
isters unchanged, load and run the Gain and Transmission Loss 
preogram according to its instructions. 


|| ۳۹ about units: 


As currently written, the programs use metric units; 
meters for depths, and m/sec for sound velocities. To con- 
E -he programs for english unit input data, feet for 
depths and ft/sec for velocities, the conversion factor 
1,852 m/nm should be changed to 6,075 f£t/nm where occurring. 
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NE wEnstructions for CZ Range Programs 


Step 


٦۲660۰66 5 


Enter 5-point SVP which 
describes the Deep Sound 
Channel. (D,=SLD, D.- 


5 
Bertem OË DSC.) 


Press R/S to perform pre- 
liminary calculations 


Enter Receiver Depth and 
Source Depth. 


Range to inner edge of 
SI REZA) and range to 
outer edge of CZ (RCZo) 
may be displayed by use 
of User Control keys B 
and C respectively. 
Ranges are in nautical 
miles. 
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Ese Allocation for CZ Range Programs 
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Step Keys Code | Explanation otep keys vode | explanation 
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Instructions for CZ Gain and Transmission Loss Program 


Step INSETUCELONS Input Keys 3116٤ 
L Run one of the CZ Range 


Programs. Leave the 
callentator on and all of 
the data storage regis- 
ters unchaged. 


2 Load the CZ Gain and 
Transmission Loss Program. 


3 Enter a zero if the Shal/ 0 A G (dB) 
Shal or Crosslayer program 
was used, or enter a one 
if the Deep/Deep program l 
was used. 


Or 


4 Enter a frequency of in- COATS QU GIBT) 
gderest ın Hz. 


5 Enter the number of a CZ n ር TEUS GB) 
annulus of interest. (Step $ 
5 may be repeated for as 
many CZ annuli as desired.) 


6 To determine the TL, 


values for different 
frequencies, return to 
step 4. 
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peerage Allocation for CZ Gain and Transmission Loss Program 
pa ters: 


2 
RO: Ar SO: 9; A: 950 / 9po 7 f 
Rl: “1 Sl: g, “ሠን ኑ "rswp 
ED: = S2: ږو‎ C: cos 9, Ve 
R3: C, S3: 94 SE 5-55 
R4: Ce S4: a terms D: cos 9 
w+1 

R5: ር S5: p ፅ 

3 E: Fo / PgR / PRR 
RO: REZT S6: 09 ٦ 
NEE, / a S7: 
R8: RCZO S8: 
R9: G terms / G S9: 


0801116-1 EE SEUSS and Use: 


mS OF, On for shallow 72۰۰0۶ ۶۰۰۱۱۱700 7 
sound Source. 


l: Off, On to decrease tne DE One, set Sy data entry 
increments used in 01 1 ٣۵ cound. 
fina ing 9.5: 


Display Status: DSP 0 


User Control Keys: 


A: Compute G a: 
B: Compute a BD: 
ea Compute TL. C: 
D: d: 
E: e 
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: Me... nn ہہ سو ہے ور رہ‎ ሮፓ 
Step Keys Code | Explanation Para, ode aict t Lor 


mg مول وسو‎ ee ee ee وس‎ መ: re ` "9 reg ` ER mn 










081 بر مل‎ cup Set FO for 851 l #1 Large 9 | 


: و‎ SLD 
00 2 ef 16-43 | Spa] S ፀ52 + ECC S 

. : USO Urea x 2 2 
882 578 [5 21 88|” 853 ro 22 55 | increment 





iTS +8 ٦ 
increment 


of Garp 














815 
816 ደአ! -21 
Bir BCL? 35 85 | 2 AT) 
818 ደርርሀ 25 14 

619  &SBa 23 16 11 
828 —ST*8 35-55 8 
821  KCLS 36 85 
822 ENT -21 
023 1 -21 





B 
E 


























824  RCL? 36 87| 2 Ar 

25 «109 2 Convert 

826  GSBa 23 16 11 Hc s 3070 

027 ETag 35-55 68 to meters 

028 RCLI 36 ði 

829 RÜLD 36 14 

938 > -24 

831 RCL: 36 45 

832 > -24 

ez Soup ' 36 14| 45 

634 SBb 22 15 12 
ı 825 X -35 | 
| 036 Stra 35-55 88 : | 

8 35-5 DS | 
| 4 STOE 5 | 
Per 106 88S  RCL2 36 82 
898 RCL 36 @1 | نټ‎ 
991 - -24 50 
Store r 098 Cis- 15 de : 
0 دجوم‎ Far 19 2304 negative if 


source shal) 
"neci ron 


Small inc. 
oí Born 

Decrease 

esLD inc. ^? 





CZ Gain and Transmission Loss Program 
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Step Keys Code | Explanation Step Keys Code | Explanation 


| ں 






















163 > -24 
184 Cds 2 
183 
Wa RCLL b 

197? . RCL4 5 84 


Second term 


181 151 4 Bs 
182 ረ 7 152 i 51 
A 


SEE 
























188 x 37 

182 ርያ ae D otore a 
11። = 

ከ5 C05 Clear S4) 
116 RCL4 6 

mie REL] o Bl 


114 = EG 
mss 


converta 
from dB/kyd 














to GB/m 
122 dog rn 
123 3 
ARCAL 
1/6 
Is 
Compute and | 178 2.6 ۶ T) 
Display G jo d 
158 N -35 
151 ን -41 
192 RCL? 36 áp |a(nRCZi) .(0 . 
133 ëlo 133 x -35 
134 EEX Convert 1 184 + -55 
وم 185 00-6“ 3 ورا‎ 36 83 | Compute & 
: 125 - -45 | Display wee 
rar se 187 3 Ei 
188 *xLBLa 21 lo 11 
Attenuation 167 STOC 53 د‎ 
Coefficient | 1% a -33 
121 (el -4 Í 
192 = “ረፋ 
Gs 193 3700 35 14 r 
First term 134 656 ረ«5 16 1ረ Sn ne 
125 RELL $6 13 
196. ፎ5ይዌይ 23 16 le 
57. - -45 
128 -39 
198 PELL Tu 
286 : -i4 





اککاقسم — A ጨሙ-ጨ.. PY‏ ےی می س و سه و وو رس ere‏ — — 


CZ Gain and Transmission Loss Frogram 
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Step Keys Code | Explanation | e E 70 











201 ٤ 36 5 | | 
83د‎ 1821 16 26 46 

283 > CIT AT Se | 

284 ጸ82 A | 

295 : gz | (Continued) | 

285 X -39 

24 ۶۸ھ 207 





Sos O Sin 


Subroutine 





- ወ ee e ፡ A — —— سا‎ 





A  መመመመ--- =>  :- 
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